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Abstract 

We formulate quenched chiral perturbation theory for heavy-light mesons 
coupled to pions, and calculate the one- loop chiral logarithmic corrections 
to /b, /bs, Bb and Bbs- We also calculate these corrections for "partially 
quenched" theories. In both theories, the chiral logarithms diverge in the 
chiral limit, indicating that (partially) quenched theories should not be used 
to study this limit. Comparing the chiral logarithms to those in QCD, we 
estimate the errors caused by (partial) quenching. By forming suitable ratios, 
we can reduce the uncertainties in our estimates. 



I. INTRODUCTION 



Numerical simulations of lattice QCD are providing increasingly accurate values for 
hadronic matrix elements. Among the most interesting is the quantity /^-Bs, knowledge of 
which allows extraction of the CKM matrix element Vtd from experimental results for B^ — B^ 
mixing. Unfortunately, present large scale simulations use the "quenched" approximation, 
in which internal quark loops are dropped, or the "partially quenched" approximation, in 
which the quarks in the internal loops have a different mass from the valence quarks. Unlike 
other sources of error, such as discretization, it is difficult to systematically reduce the errors 
due to quenching. Thus it is important to estimate these errors, and in this paper we do so 
for the decay constants and B-parameters of heavy-light mesons. 

We estimate quenching errors by comparing the size of the chiral logarithmic contribu- 
tions to the quantities of interest in QCD and quenched QCD (QQCD). Similarly, errors 
due to partial quenching are estimated by comparing chiral logarithms in QCD and in par- 
tially quenched QCD (PQQCD). Chiral logarithms arise from the long-distance part of the 
clouds of pseudo-Goldstone bosons (PGB's) surrounding all hadrons. These clouds are rad- 
ically different in the two theories, there being, for example, a long-distance rj' cloud around 
quenched hadrons, while there is no such cloud around hadrons in QCD. Thus, in gen- 
eral, there are significant differences in the chiral logarithms. In particular, we find that fs 
and Bb calculated in QQCD should diverge in the limit that the light quark mass vanishes, 
whereas in QCD they should approach a finite limit. Such divergences are also predicted for 
light meson quantities, e.g. /^r, and ml /nig 0,0, and, indeed, have been observed 0. 

Our estimates of quenching errors in individual quantities, such as Jb, are necessarily 
crude, because chiral perturbation theory tells us nothing about analytic contributions. We 
can do much better, however, by considering ratios of quantities designed so that all unknown 
contributions vanish. The only corrections to these ratios come from 0{mg) terms in chiral 
perturbation theory. 

Our calculation is done by combining two technologies: the chiral Lagrangian for heavy- 
light mesons ^ (which incorporates the heavy quark symmetry), and chiral perturbation 
theory for the quenched and partially quenched approximations [0,0.0 Our results are thus 
valid up to corrections of 0{1/Mb), which we expect to be small. 

The outline of this paper is as follows. In the following section, we explain how to extend 



quenched chiral perturbation theory to include heavy-light mesons. In section |T| we outline 



the calculation and give our results. We discuss our estimates of the quenching errors in sec. 



rVl , and offer some conclusions in the final section. 

After our calculations were completed, a closely related paper by Booth appeared Q . He 
calculates the chiral logarithmic corrections to fs and Bb in QQCD, as well as the corrections 
to B masses and the Isgur-Wise function. Our results for QQCD agree with his, except that 
we have included two coupling constants that he overlooked. This difference is not crucial, 
however, as these couplings are expected to be small — indeed, in our numerical results we 
assume that they vanish. We keep the discussion of technical details to a minimum, and 
stress results that are new. In particular, we have (1) devised ratios in which unknown 
corrections of 0{mq) cancel, allowing a more reliable estimate of quenching errors (this 



-'^For recent reviews of quenched chiral perturbation theory and its apphcations, see Refs. 
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requires us to extend the calculations in QCD to the case of non-degenerate quarks); (2) 
considered wider ranges of the poorly known parameters in the chiral Lagrangian; and (3) 
extended the estimates to partially quenched QCD. 

In a second work, Booth has extended his calculations to include 1/Mb corrections [p!0[ . 



This introduces a number of additional unknown constants, making estimates of quenching 
errors more difficult. It appears likely, however, that the 1/Mb corrections will not be 
large enough to nullify the conclusions drawn from the infinite mass limit assumed here. A 



detailed discussion of these corrections in QCD has been given by Boyd and Grinstein [ITT 



II. FORMULATION OF QUENCHED CHIRAL LAGRANGIAN INCLUDING 

HEAVY-LIGHT MESONS 

We begin with a brief review of the quenched chiral Lagrangian introduced in Ref. 0. 
Following Morel |T2|, a functional integral formulation of QQCD is obtained by adding, for 
each light quark (a = u,d, s), a. ghost (bosonic) quark with the same mass. The fermion 
determinant of the quark sector is then canceled by that of the ghost sector. In the limit of 
vanishing quark masses, and ignoring the anomaly for the moment, the QQCD Lagrangian is 
invariant under the graded group f/(3|3)i x ?7(3|3)/j, which mixes quarks and ghost-quarks. 
Numerical evidence strongly suggests that this symmetry is broken dynamically down to 
its vector subgroup, just as in QCD. This gives rise to 36 Goldstone mesons, 9 each with 
compositions (pab ~ ga^fe, 4>ab ~ qaQb^ xib ~ laQb and Xab ~ qaqb- The first 18 are bosonic, 
the others fermionic. 

The next step is to construct an effective Lagrangian describing the low energy {p <C 
~ 1 GeV) interactions of these particles. We first collect the Goldstone fields into an 
hermitian 6x6 matrix. 



(2.1) 



which transforms linearly under the vector t/(3|3), but non-linearly under the full chiral 
group. Its exponential, however, transforms linearly 

E = exp(2i$//) , S — > L^R^ , (2.2) 

where L G U{3\3)l and R G U{3\3)r. The lowest order Lagrangian respecting the U{3\3)l x 
f/(3|3)/j symmetry is 

A™ = 4rStr(5^S9^St) + Ao str(MS + MS^), (2.3) 



where the mass term is 



and the supertrace str appears because of the graded structure of the symmetry group 
/ and Ao are undetermined bare constants. In our numerical estimates, we set the scale in 



M = dia:g{mu,md,ms,mu,md,ms) , (2.4) 
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QQCD (and PQQCD) by taking /Q^^d _ jpqqcd _ jQCb ^ ^gg MeV. (This is shghtly 
less than = 132 MeV because of the extrapolation to the chiral limit.) Having made this 
choice, we use the common symbol / for all three quantities. In the same spirit, since Aq 
does not appear in the final expressions, we do not distinguish between its versions in the 
three theories. 

The next step is to account for the anomaly. The chiral symmetry of the quantum theory 
is, in fact, [SU{3\3)l x SU{3\3)r]xU{1). In particular the field 

(fj' — f]') str$ 

$0 = ' = oc str In S = In sdetS (2.5) 

V 2 v6 

[rj' is the normal SU{3) singlet meson, rj' its counterpart in the ghost sector) is invariant 
under this smaller group, but not under the full U{3\3)l x 6^^(313)^. Thus we can multiply 
each term in the 0{p'^) quenched chiral Lagrangian by an arbitrary function of $o- The 
full structure is discussed in Ref. 0; in our one-loop calculations we need only the terms 
quadratic in 

Cnght = Unv + «(5;.$o)' - val^l ■ (2.6) 
We are using the notation and conventions of Ref. [Q (in particular f^^ = 132 MeV), except 



that we use mo instead of fi, following Ref. [|I4 



It would appear that a and mo can be absorbed into the quadratic part of Cmv, renor- 
malizing the fields and masses of the 1]' and rj' particles. In fact, as explained in Refs. 0,0, 
a and mo should be treated as interactions. Diagrams with multiple insertions of either 
coupling on an rj' or rj' line cancel (because they correspond in the underlying theory to 
diagrams containing quark loops). Thus the mo term does not lead to a mass for the rj' or 
rj'. They remain light, and must be included in the effective Lagrangian. This is in contrast 
to QCD, in which the rj' gets a mass of 0(mo) and can be integrated out. In our calculation 
(as in all work in quenched chiral perturbation theory to date) we assume that a and ml 
are small, so that we can work to one-loop order. 

Diagonalizing the mass term in Cinv, one finds that the flavor diagonal mesons have 
compositions uu, dd and ss, with masses 



2 8Aom„ 2 SAomrf ^ SX^rUs 



respectively. The flavor non-diagonal particles have the same masses as in QCD, e.g. m|-+ = 
4Ao(mtt+ms) / p. The states containing ghosts have the same flavor compositions and masses 
as those without. 

This ends our review of the construction of Bernard and Golterman. We now extend the 
effective quenched Lagrangian to include the ground state mesons containing a heavy and 
a light quark. As in the discussion above, the result is a simple extension of that in QCD, 
which was worked out in Refs. . The only differences are the presence of fields including 
ghost-quarks, and extra interactions involving $o- 

In the heavy quark symmetry limit the pseudoscalar and vector heavy-light mesons are 
degenerate. They are destroyed, respectively, by the fields and Pa^^'' ■ These fields have 
an implicit label v denoting the four-velocity of the corresponding particles. In terms of the 
underlying quarks the flavor composition is -P^*^-* ~ Q^a^ so that the meson destroyed by 
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p{Q) contains a heavy quark Q and a light anti-quark of flavor a. In QQCD this index runs 

over u, d, s and u, J, s. Thus -P^*^^s are commuting fields, while are anticommuting. 

In the chiral limit all flavors are degenerate. 

To display the heavy quark symmetry we form the spinor field 



1+ 



p*{Q)^tJ- _ piQ) 



75 



which transforms as 



(2.8) 



(2.9) 



where S is an element of the heavy quark spin SU{2)y group ||T5[. We do not elaborate on 
this symmetry here, since the part of our construction which uses this symmetry is identical 
to that for QCD. To display the transformation properties under SU{3\3)l x SU{3\3)r we 
follow the standard procedure of introducing U G SU {3\3), defined such that 



e = exp(z$//) LiU^ = U^R^ , 
where now L G SU{3\3)l and R G SU{3\3)ii. Then H transforms as 

a 

The conjugate field is defined by 

riQ) 
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H 



iQ) 



{UabHf^)S-' 



(2.10) 



(2.11) 



(2.12) 



This field creates heavy-light mesons of velocity and has nothing to do with the anti- 
particles of these mesons. 

Now we can write down the part of the effective Lagrangian containing Hf^^ 



ba 



ba 



+ iz g/ straTvHf^H^^h,!^ str S}d'^ - ^d^( 



+ 



(2.13) 



where the ellipsis denotes terms of higher order in momentum, quark masses, and l/mq. 
stfa indicates that the supertrace is to be taken over the repeated index a. Summation over 
other repeated indices is assumed. Finally, the trace Tr is taken over the Dirac indices. Cq 
is invariant under the chiral group, the heavy quark spin symmetry, and, of course, Lorentz 
transformations. In principle, each term in Cq can be multiplied by a real and even but 
otherwise arbitrary function of $o-0 We do not display these functions as they contribute 
only at two-loop order. Aside from these (implicit) functions, the only difference between 
the quenched Lagrangian Cq and that in QCD is the appearance of the g/ term. Since 



^Note that a possible term linear in $0 of the form ^oStraTrHaHaj5 vanishes identically. 
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str ^^d^C, — ^d'^^^ oc 9^$0) this term gives rise to a coupling of the t]' and t]' to the heavy- 
hght mesons. Though such a couphng exists in QCD, it is not needed in the QCD chiral 
Lagrangian because the rj' is integrated out. 

To discuss i?-parameters, we also need the Lagrangian describing the interactions of 
PGB's with the antiparticles of the heavy-light mesons, i.e. those with composition QaQ- 
We follow Ref. [0 and introduce a field which destroys such mesons 



This transforms under heavy quark spin symmetry and SU{3\3)l x 5't/(3|3)/j as 

Hp iU.,Hp)S-' . 
The conjugate field H^'^^ is defined as for H^'^K The invariant Lagrangian is 



(2.14) 



(2.15) 



Cq = -i stfaTrv ■ dHpHp - \istraTrHpHf\^' [s}d^'i + ^9^^^ 
+ hgfstraTrHpHf\,^,str l^d^^ - ^^'^^^1 + ■ • • • 



ba 



(2.16) 



The total quenched chiral Lagrangian is then C = Cught + + C.q- 

We are interested in decay constants and 5-parameters. These are defined, respectively, 

by 



(0|u7^(l-75)6|5 {v)) = -ifBP^, = -ifsmBV^ 



(2.17) 
(2.18) 



with similar equations for mesons of other flavors. We thus need to know how the left- 
handed currents, = qal^iX~l^)Qi ^-^d the four-fermion operators composed of two such 
currents, can be written in terms of the fields in the quenched chiral Lagrangian. The result 
is that the matrix elements of the current involving the Q field are given by 



(2.19) 



where the ellipsis denotes non-leading operators in the chiral and heavy quark expansions. 
The "potential" Vl can contain both odd and even powers of $0) with arbitrary coefficients — 
and we find that the linear term does contribute at one-loop order. We normalize the 
potential so that Vl(0) = 1.0 Then, taking into account the factor of ^niQ-q absorbed into 



•^We note in passing that similar potentials occur in the chiral representation of currents composed of light 
quarks. They do not, however, give rise to one loop corrections to and Jk- This is because CP, together 
with quark-interchange symmetries, require that the potentials be even functions of $0. 
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the fields P^'^\ we have k = foq^/rn^ at leading order. All this is just as in QCD p, except 
that the index a now runs from 1—6. 

Similarly, following Ref. |]16|, the mixing matrix elements of the AQ = 2 operator Oaa = 
Qjfi{l — 75)gaQ7'^(l — 75)^0 (no sum on a) can be represented in the effective theory by the 
operator 



Oaa = Y.(^TrUHf\,{l - 75) TrUHi'^Vil - l5)Vo{%) + ■ ■ 

b,c 



(2.20) 



Again, the arbitrary potential is normalized so that Vo(0) = 1, and can contain terms linear 
in $0 which will contribute at one-loop order. This representation is not unique, but, as 
explained in Ref. [|TB[, other choices are equivalent since they all lead, after evaluating the 
Dirac traces, to the form 



(2.21) 



At leading order /3 = ^f^mBBs- 



A. Partially quenched theories 

Partially quenched theories are a first step towards removing the quenched approxima- 
tion. In these theories, not all the quarks are quenched, so that there are more quarks than 
ghost quarks. The particular example we focus on is that used in recent simulations (e.g. 
see Ref. fl^]): the lattices are generated with Nf quarks of mass ruf, while valence quarks 



of mass rriq are used to calculate fs^ and Bb^- Although Nf = 2 in present simulations, 
we keep Nf as a free parameter, so as to make our results more general. In the limit that 
rriq = mf the theory is unquenched and QCD-like, having Nf flavors of degenerate quarks. 
Thus we refer to it as partially quenched QCD. 

It is straightforward to generalize the development of the previous subsection to PQQCD. 
Most of the work has been done by Bernard and Golterman, who developed chiral pertur- 
bation theory for partially quenched theories 0. There are only two alterations. 

• The index a labeling quarks and ghost-quarks now takes on Nf + 2 values: a = 
(g, /i, . . . , /jvy., q). (We could introduce quenched u, d and s quarks, as in our discussion 
of QQCD, but this does not affect the results for the quantities of interest.) The 
initial symmetry is thus U{Nf + 1|1)l x U{Nf + 1|1)_r, the fields S and ^ being 
(A'^^ + 2) X [Nf + 2) matrices. The form of the heavy-light meson Lagrangian is 
unchanged. 

• The a and rriQ vertices can be iterated in the SU{Nf) sector, so that the mass of the 
1]' in this sector is shifted from the quark mass contribution, m^j = SXofnf/ f"^, to 



ml, = —r^ ■ 2.22 



The qq meson mixes with this rj' through the and a vertices, leading to a propagator 
with a comphcated pole structure, given explicitly in Ref. 0. 



7 



/Hi 





(a) (b) (c) 





(d) (e) (f) 

FIG. 1. Graphs contributing to the one- loop renormalization of the decay constant. Thick and thin 
hnes are propagators for heavy-light mesons and PGB's respectively. Solid circles are vertices proportional 
to g or g/, crosses are insertions of the interactions proportional to a or mo, and the filled squares denote 
the left-handed current. 

To simplify the resulting expressions, we assume that both mf and are small, and that 
Nj is large enough that the SU (Nf) rj' is heavy. We then have a hierarchy of scales 

~A, (2.23) 

and we can drop terms suppressed by rn?^^/ {miNf /?>) or m'jj/^mlNf /3). When this is done, 
all contributions that depend on the "quenched couplings" (i.e. mg, a, gf, V^(0) and Vq^O)) 
turn out to be higher order in the chiral expansion than those we keep. It should be borne 
in mind, however, that, for the largest values of mjj: used in present simulations, corrections 
proportional to mjf/ml may be substantial. 



III. CALCULATION OF CHIRAL LOGARITHMS 

We first calculate the one-loop chiral logarithms contributing to /g^. Strictly speak- 
ing, we calculate the logarithms for the quantity Ksy/mq-s = /gsmg^. At leading order in 
l/mg, however, these results apply also for the decay constants, because chiral logarithms 
in mQs are suppressed by l/mq. The graphs which contribute are shown in Fig. (|l]). The 
calculations are standard — the details are as in QCD [jl6l. In particular, the heavy-quark 



symmetry allows us to choose the external meson to be either the P or P*. We choose the 
former, which allows us to extract the result by contracting the external current with v^. 
The three-point vertices in Cq connect a P only to a P*, so that the heavy-light meson in 
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the loops is always the P*. This means that graphs (a) and (b) vanish identically (as in 
QCD), because the factor of from the current annihilates the P* propagator. 

The absence of quark loops in QQCD manifests itself as cancelations between inter- 
mediate states containing ghosts and those without. Thus, while graph (c) can give rise to 
contributions proportional to g"^, g/^, and gg/, only the latter is non-vanishing. Similarly, the 
only non-canceling contribution from (d) is that proportional to g"^, and the contributions 
from (e) cancel entirely. The net effect is that only loops involving the ss PGB contribute — 
just as expected since the s quark is the only light quark available. Expressing the decay 
constant as 

^QQCD ^ 

(/q.)^^^° = ^ 1 + (A/q.)QQ^° + c?{A)ml + c^iA)ms + • • • , (3.1) 
we find the chiral logarithm to be 



\QQCD 



g (1 + ^9q) (2am3, - rngj - I ^ 



log 



mi 



A2 

(3.2) 



Here m^^ is defined in Eq. ( ^.7[ ), and A is the ultraviolet cut-off on the logarithmic diver- 
gences. The A dependence in A/ is canceled by that of the analytic terms proportional to 
c^2 in Eq- |3-1| - These constants result from higher order terms in the expression for the 
current L'^ and in the heavy-light meson Lagrangian. Their structure is discussed by Booth 

i- 

The quenched result for fq^ (g = u, d) is obtained from Eq. ^]2] simply by substituting 
iriqq for iJiss- Present lattice simulations, however, do not directly work with physical u and 
d quarks, but instead extrapolate from heavier quark masses. We mimic this procedure 
analytically in the next section. 

The entire result in Eq. |3.2| is a quenched artifact, being proportional to one of a, mg, gl 
or V^(0). The quenched logarithms are thus very different from those in QCD. In particular, 
the term proportional to diverges in the chiral limit, whereas the chiral logarithms in 
QCD, to which we now turn, vanish in this limit. 

The chiral logarithms have been calculated in QCD in the limit iriu = md = in Ref. 
||16|| . We need, however, the result for the general case of three non-degenerate quarks. We 
use a similar expansion to that for QQCD 

^QCD 

(/q^)'^''° = ^= 1 + (A/q^)Q^° + Ci(A)(m„ + m, + m,) + C2(A)m, + . . . , (3.3) 



where q = u, d or s. Note that the analytic terms have a somewhat different form to those 
in QQCD. The chiral logarithm comes from Figs. |1] (c) and (e), and is 



2 



(A/q,)^^° = - (^(1 + 3/) T.iT!),,ml log ^ . (3.4) 

The sum runs over the PGB's (eight in QCD) in the mass eigenbasis, with being the 
corresponding group generator, normalized such that tr[TaT{,] = ^6ab- This formula hides 
the fact that, for 7^ m^, the and rj mix, so that one must diagonalize the mixing 
matrix to determine the appropriate T^'s to use. 
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(a) (b) 




(c) (d) 
FIG. 2. Graphs contributing to the one- loop renormahzation of Bg-q. The sohd square denotes the 
operator Oaa- Other notation as in Fig. 1. 

The result for PQQCD is obtained from the same graphs as that in QQCD, except that 
the crosses in Fig. [l|(b) and (d) now include a projection against the SU{Nf) rj'. The 
expansion for the decay constant is 



^PQQCD 
_nPQQCD _ ^ 



1 + (A/q,)™° + cf«(A)m; + c^^{A)m, + ...], (3.5) 



with the chiral logarithm being 

(Afej-^^" = -^(1 + 3gl^) m^, log ^ + ^(mj, - 2mJ,) log ^) . (3.6) 

Here rugf is the mass of the meson composed of a valence quark (mass rUg) and a dynamical 
antiquark (mass mf), and to the order we work is given by rrigj = {rrigg + m'jj)/2. A check 
of Eq. pj.6| is that, in the limit rrig = mj, it equals that for QCD with Nj degenerate flavors 
(easily obtained from Eq. |3.4| ). We stress again that terms proportional to a, and gi are 
absent. 

We now turn to results for i?-parameters, the relevant graphs for which are shown in 
Fig. 1^. We define Ai? in the same way as A/ (Eqs. P?T| , |^ and ^.5] ) — one can show that 
the analytic terms have the same form as for decay constants for all three theories. The 
values of the constants for 5-parameters are, of course, not related to those for decay 
constants. Our results for the chiral logarithms are 



\QQCD 



(47r/)^ 



2 , / 2 2a o mi 



' log^, (3.7) 
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(4vr/)2 
1 



;i-3/)E2[(Tj,/m^log^ 



pq) 



1 



log 



m 



99 



A2 



(3.8) 
(3.9) 



Only flavor-neutral PGB's contribute in all three theories. This is explicit for QQCD and 
PQQCD, while for QCD it follows because the expression picks out generators whose gg-th 
element is non-vanishing. As required, the PQQCD result reduces to that for QCD with Nf 
degenerate quarks when mf = niq. 

The quenched result for AB is similar to that for A/ (Eq. p.2|) , except that the first 
term is not a quenched artifact. It is of the form rn^glogrrigq, but is not proportional to a. 
Indeed, this term is exactly the contribution that would be obtained from the QCD formula, 
Eq. p78| , if one were to use the quenched spectrum of flavor-neutral mesons, i.e. states with 
composition uu, dd, and ss. In other words, the quenched result consists of that from QCD, 
modified because the spectrum of PGB's is different, together with contributions which are 
quenched artifacts and involve the rj'. This is the same situation as for Bk, and can be 
understood by an analogous argument using quark-line diagrams 



IV. ANALYSIS AND DISCUSSION 

To use our results we need to know the numerical values of the dimensionless parameters 
mo//, a, Qq, gpQ, gf, fVl{0) and /V^(0). These are parameters of the (partially) quenched 
theory, and can, in principle, be determined by simulations. At present, however, these 
parameters are not well known, so we allow them to vary in "reasonable" ranges, determined 
as explained in the following. To compare to QCD we also need to know g. 

There is now evidence for a non-zero value of rrio/f The result is conventionally 

expressed in terms of 6 = 2ml/ (3(47r/)2). A reasonable range for 6 appears to be 0.05 — 0.15. 

The evidence for a is less clear. Like mg, it is proportional to 1/A^c, and thus, hopefully, 
small. The simplest choice is to assume a = 0,^. A better estimate, however, can be 



obtained from the data of Ref. [113] on the r]' two-point function. Fitting the data to the form 
expected from Cught (Eq- |2.6D,one finds a ~ 0.7 and uiq ~ 720MeV {6 = 0.13). This result 
may have large discretization errors since it is obtained at a relatively large lattice spacing 
(/3 = 5.7). Nevertheless, it suggests a positive value for a. We assume that a lies in the 
range — 0.7. 

There are no lattice results concerning gq or gpQ, so we assume that the value is similar 
to g. Experimental results for r{D* —>■ Dtt) imply 0.1 < g^ < 0.5 (see the discussion in Ref. 
0). We use this as our reasonable range, and take gq = gpQ = g"^. To be consistent with 
Ref. p we use g^ = 0.4 as our preferred value. 

Finally, we come to g/, V/,(0) and Vq{0), about which even less is known, g/ is the addi- 
tional p('3)-P^('3)-ry/ coupling arising from the fact that the qq pair in the T]f can annihilate 
into gluons. It follows that g'/gq oc 1/Nc, though we have no information concerning the rel- 
ative sign of gq and g/. We consider two possibilities: a conservative choice, —gq < g/ < gq, 
and the more likely situation in which g/ is small, \gf/gq\ <^ 1. Similarly, V^(0) and Vq{0) 
correspond to additional couplings of the i]' to the current and four-fermion operators, re- 
spectively, involving annihilation into gluons. Thus both are of 0{1/Nc). In the following. 
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FIG. 3. Quenched chiral logarithmic contribution to fg-^. 

we assume that both vanish. In fact, g/ and V^(0) appear in the same way in both the 
correction for the decay constant and 5-parameter, so that V^(0) can be absorbed into g/. 

With these parameters in hand, we begin by plotting [Afg-s)^^'^^ (Eq. |^]^), as a function 
of m^g. Here "s" indicates a generic quenched "hght" quark. For the couphngs we take 
6 = 0.15, gq = 0.4, A = 1 GeV, and display the results for a = 0,0.7 and gf = —gQ,0,gQ. 
The following important features turn out to hold for all values of the couplings in our 
reasonable ranges. 

• The correction diverges as rriss —>■ 0. This reflects the unphysical cloud of light r]' 
particles surrounding the quenched heavy-light mesons. Observing this divergence in 
numerical results would give further confirmation of the validity of quenched chiral 
perturbation theory. In practice, this may be difficult, since the divergence does not 
set in until rriss falls below ~ 100 MeV. 

• For masses heavy enough to be away from the divergence, the correction can be sub- 
stantial. This is the region in which most present simulations are carried out: "light" 
quarks range in mass from roughly to mf^^^, so that 0.17(GeV)^ < m^^ < 
0.5(GeV)^. For these masses, and for couplings in our reasonable ranges, we find that 
the magnitude of the correction can be as large as 30%. If \g'/g\ ^ 1, however, the 
correction is smaller, not exceeding 15%. 

• The departures from linearity in the range 0.17(GeV)^ < m^^ < 0.5(GeV)^ are fairly 
small, particularly for g/ > 0. In the following, we mimic numerical simulations, and 
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FIG. 4. Chiral logarithms for Nf = 2 dynamical quarks: PQQCD with to/ fixed at two values, and the 
QCD-like theory with to/ — TOs. Parameters are A = 1 GeV and gpq = 0.63. 



define a quenched "/b" by linear extrapolation from the masses m^^ = 0.25 (GeV)^ 
and 0.5 (GeV)2. 

• Afq-g decreases with increasing m^g. Using the definition of fs just given, the chiral 
logarithm alone thus implies fBs/fs < 1 in QQCD. This is in contrast to QCD, where, 
for g = 0.63 and A = 1 GeV, (A/^JQ^d = o.38 > (A/b)Q^° = 0.21. We return to 
this point below. 

In Fig. ^ we show the size of the chiral logarithms in partially quenched QCD. We 
have chosen two values for the dynamical pion mass which are in the range used in present 
simulations, mjj = 0.15 and 0.3GeV^. The divergence at small m^^ remains, although now 
it is proportional to mjj: rather than ttt-q. Again the corrections are substantial, although 
the dependence on m^^ is weaker than in QQCD. We also show results for the unquenched 
two flavor QCD-like theory obtained by setting rrif = rris- Although there is no divergence 
at small m^s, the curve is not very different from those for PQQCD in the mass range of 
typical simulations. 

An example of the results for (ABqs)^^'^^ , is shown in Fig. ^. There is very little 
dependence on a, so we show results only for a = 0.7. Similar comments apply as for A/q^, 
except that 

• The sign of the term proportional to g/ has changed; 

• The divergence at small m^^ is weaker. This is because it is proportional to 1 — Sgq, 
to be compared to 1 + Sgq for Afq-g. The sign of the divergence depends on gq. 
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FIG. 5. Quenched chiral logarithmic contribution to Bq-^. Parameters as in Fig. 3. 



A. Quantitative estimates of quenching errors 

These results suggest that quenched artifacts give important contributions to and 
Bb- They do not, however, allow a quantitative estimate of quenching errors. For this we 
need to know something about the non-logarithmic contributions. For example, to estimate 
jQQCD^jQCD clearly need to know both and the difference between the 

analytic 0{mq) and ©(mg) corrections. These can compensate for differences in the chiral 
logarithms. 

The cleanest way of proceeding is to consider ratios which cancel as many of the unknowns 
as possible |jl9|. Single ratios such as fss/ Ib (with f^^'^^ defined by extrapolation as 



explained above) do not depend on /-QQCD^^qcd^ depend on the difference of 0{mq) 

terms in QCD and QQCD. To remove this dependence, we must form a double ratio, which 
requires the use of three non-degenerate valence quarks. Calling these u, d and s, and 
assuming m„ < rrid < rUs, we construct 

g^^(hE]'^(fSl]'^- .1 (4.1) 

2(,n,-,n.)— ^fe) + 0(m?). (4.2) 

and the similar ratio Qb involving the 5-parameters. To the order that we work in chiral 
perturbation theory, the quark masses appearing in Q j can be replaced by the corresponding 
squared meson-masses, e.g. m^^. Any dependence of the form A/g^ oc a + &(m„ + 
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TJid + Tils) + cniq cancels in Qj, up to corrections of O(m^). In particular Q/ is independent 
of the cut-off A. Note also that Qf has a well defined limit as nid nig or — m„. An 
interesting special case is if the masses are in the ratio m„ : irid : = 1:2:3, for which 

Qf = foufos/Uodf'- 

Qf and Qb are thus quantities which we can calculate with some reliability in all three 
theories which we consider. The differences give a direct measure of the (partial) quenching 
error. The results for QQCD follow from Eqs. |3.2| and |3.7| 



Q 



QQCD 
/ 



^QQCD 



(4vr/)2 
(4vr/)2 



\ V6 



(4.3) 



2a 



l-34)(— -1) 



Ri-5 



(1-3^7^) 



R2 



(4.4) 



where Ri^2 are functions of only the quark mass ratios 

mslog{ms/md) m„log(md/m„) 



Ri 
R2 



2ms 
rris - m,, 



TUs -TUd TUd- rriu 

\og{ms/md) log(mrf/m„) 



frid 



TUd - TUu 



(4.5) 
(4.6) 



For QCD, after some algebra, we find 



^QCD 
^QCD 



:i + 3^7^) rris 



m,, 



The vanishing of Q^^^ seems to be accidental — we expect the contributions of O(m^) to be 
non-vanishing. Finally, for PQQCD, we have 



ml+ log(m|.+ / ml+ ) mj^o log(m^o ) 



rris - nid 



(4.7) 
(4.8) 



Q 



PQQCD 
/ 



:i+3^?i>q) 



Q 



PQQCD 
B 



(47r/)2 

;i-3(? 



N. 



1 



-^m^^R^ + ■^^^mjj-R2 - —m^^Ri 



pq) 



(47r/)2 

where the additional ratio is 



2 



2N 



f 



N 



f 



R, 



4m o 



rris - rud 



md + rrifj rrid - rriu 



(4.9) 
(4.10) 

(4.11) 



To make numerical comparisons we choose two sets of quark masses: m„ : md : = 
1:2:3, ratios that are close to those available from a typical simulation; and : rrid : 



1 : 2 : 10, ratios now attainable in QQCD for 



phys 



IQ, and which approach 



those of the real world. The results are shown in Table |. Typical simulations have x ~ 1, 
and we expect y to lie in the range 1/3 — 1. For the smaller mass ratios, we see that there 
is a difference between in QCD and QQCD, but the quenching error is likely not larger 
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TABLE 1. Results for ratios Q in QCD and QQCD, in terms of a; = m^jO.SGeV^ and y = (5/0.15, with 
g/^Vl{0) = VSiO) = 0. 



THu : 


md : ms 


9'=9h 


^QCD 


gQQCD 


^QQCD 


1 : 2 


: 3 


0.4 


0.009 X 


0.024a x + 0.024y 


0.0042; (1.5-a) ■ 


1 : 2 


: 10 


0.4 


0.026 X 


0.082ax + 0.182y 


0.015x(1.5-a) ■ 


1 : 2 


: 10 


0.16 


0.017X 


0.055ax + 0.123y 


-0.039x(1.5-a) 



TABLE IL Results for ratios Q in PQQCD, with A^/ = 2 and gpq = 0.63. 



niu ■■ 


md ■■ ms 


mjf = O.lSGeV^ 


= 0.5GeV2) 
mjf = 0.3GeV2 


mjf = O.lSGeV^ 


= 0.5GeV2) 

mjf = 0.3GeV2 


1 : 2 
1 : 2 


: 3 
: 10 


0.021 
0.129 


0.035 
0.216 


-0.0015 
-0.012 


-0.003 
-0.024 



than a few percent. This is probably comparable to the statistical errors in the best present 
simulations. The corresponding error for is much smaller, and unlikely to be measurable. 
Other values of the couplings lead to similar conclusions. 

For the larger mass ratios, there is a substantial quenching error, possibly as large as 
20% in Qf. This is a reflection of the divergence at small rriss in QQCD. The quenching 
error in is now larger, but, as shown by the last column, is quite sensitive to qq. 

We view these estimates as lower bounds on the quenching errors in decay constants and 
i?-parameters themselves, since the latter errors might largely cancel in the ratios defining 
Qf^B- The estimates imply that one should not simulate QQCD with quark masses that are 
too small. They also suggest that quenching errors in i?-parameters are smaller than those 
in decay constants. 

Similar conclusions hold also for "partial quenching" errors. The results for qPQQCD 

are 

not linear functions of m^^ and thus more difficult to display. The values change little in 
the range m^^ = 0.25 — 0.5(GeV)^, however, so we quote the value for m^^ = 0.5(GeV)^. 
There are no such problems for Q^^'^^ , since, for Nj = 2, it is independent of m^^ for fixed 
quark-mass ratios. Our results are given in Table ||. They are much closer to the QQCD 
results than those from QCD. For these quantities, partial quenching appears to be of little 
help. The only caveat is that, as described in Sec. |11 A| , we have neglected terms suppressed 
by {3/Nf){m'jj/ml), which may not be small for the larger values of fn'j^. 

It would be interesting to test the results for Q in quenched simulations. One way of 
doing this would be to plot Qf/R2 versus m^gi?i/i?2. The results for all mass ratios should 
lie on a universal straight line, the intercept of which is proportional to 6. An important issue 
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facing such a test is the size of higher order chiral corrections. To address this, we assume the 
correction effectively adds to A/ a contribution of size (mss/lGeV)^. For mj.^ = O.SGeV^, 
this is a 25% correction to the decay constant itself. The contribution to Qf turns out to 
be independent of 

The importance of this correction depends on the quark masses. Taking m^^ = O.SGeV^, 
one has 5Qf = 0.06 for m„ : : = 1:2:3, which is larger than the supposedly leading 
terms. On the other hand, 6Qf = 0.10 for m„ : : = 1 : 2 : 10, which is smaller 
than the leading terms. It seems likely, though, that one would need to work at somewhat 
smaller quark masses to suppress such 0{m'^) terms. 



B. Qualitative estimates of quenching errors 

We turn now to the quenching errors in the individual quantities fs, /b^; and Bb^- 
The notation we use is exemplified by 

A(A/5) = /QQ^°//Q^°-l. (4.13) 

In the previous subsection, we used QCD to refer to all theories with three light quarks, 
irrespective of their masses. In this subsection, however, QCD means the theory with 
physical quark masses. To make this distinction clear, we use fs^ instead of fg-g (and Bb^ in 
place of Bqs), although the two quantities are the same. We define the (partially) quenched 
Jb and by linear extrapolation in the light quark mass, as explained above. 

We estimate quenching errors by assuming that the dominant difference between quan- 
tities in QQCD and QCD comes from the chiral logarithms. In other words, we assume, for 
some reasonable choice of A, that 

AiAfs) ^ (A/b)Q«^°(A) - (A/b)Q^°(A) , (4.14) 

and similarly for the other quantities. This assumption is likely to be wrong for individ- 
ual quantities, because differences in chiral logarithms can be compensated by differences 
in analytic contributions. (By analytic contributions we mean not only the 0{mq) terms 
proportional to Ci^2, but also the lowest order terms, e.g. K^^'^^^.) One might expect there 
to be some compensation because the two theories are matched by equating, say, f-„, and 
so the difference in chiral logarithms for this quantity are, perforce, compensated. But it 
seems unlikely to us that this compensation will occur uniformly; for one thing, the chiral 
logarithms due to t]' loops are absent for some quantities (e.g. /^), while present for oth- 
ers. Thus we view our estimate of the quenching error as a rough upper bound (because 
there may be some compensation), which is likely to work better on the average over several 
quantities. 

We have examined the resulting estimates for our reasonable ranges of parameters, and 
for the choices A = rrip and A = 1 GeV. To illustrate their size we give in Table fT| the 
numerical expressions for the chiral logarithms in terms of a, 6 and g/, fixing g = Qq = 0.63. 
In Figs. ^ and |^ we plot the quenching error versus gf for two sets of parameters chosen to 
indicate the range of variation in the estimates. 

We draw the following conclusions. 
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TABLE III. Numerical results for chiral logarithms A*^'^^^ and A*^^^, the latter being in parentheses. 
Our estimate of the quenching error is the difference of these two expressions. We assume V^(0) = Vq{0) = 0. 



Qty. A = 1 GeV 

fs, -0.0950 + 0.385 + 0.255/ (0.38) 

fs -0.0950 + 1.145 + 0.255/ (0.21) 

fsJfB -0.00-0.765 + 0.05/ (0.17) 

Bb, -0.025 + 0.0170- 0.0665-0.495/ (-0.017) 

Bb -0.025 + 0.0170 - 0.205 - 0.495/ (-0.007) 

Bbs/Bb -0.0 + O.O0 + 0.135 -O.O5/ (-0.010) 



A = 0.77 GeV 
-0.0230 + 0.0935 + 0.0605/ (0.23) 
-0.0950 + 0.855 + 0.255/ (0.14) 
-0.072O- 0.765 -0.195/ (0.09) 
-0.006 + 0.004O - 0.0165 - O.I25/ (-0.010) 
-0.025 + 0.017O - 0.155 - O.495/ (-0.005) 
0.019 - 0.013O + 0.135 + O.375/ (-0.005) 



• There are no values of the parameters for which the quenching errors are small simul- 
taneously for all four quantities, where small means |A(AX)| < 0.1. Typically, two 
or more of the estimates are larger than 0.2 in magnitude, some being as large as 0.4. 
Viewing our estimates as upper bounds, we conclude that it is possible that there are 
substantial quenching errors in one or more of the quantities. 

• It is interesting to consider the ratios /bJ/b and BbJBb- We expect our estimate 
of the quenching errors to be more reliable here, since there is no dependence on the 
leading chiral parameters, e.g. /^QQCD^^qcd might also hope that quenching 
errors are smaller. Indeed, there are cancellations in the chiral logarithms, as one can 
see from Table pT[ (The exact cancelation of the coefficients of a and gi for these 
ratios when A = 1 GeV is not a general feature — it is a chance result, valid only for 
the particular masses we use to extrapolate to the S-meson.) Nevertheless, as is clear 
from Figs. |^ and |^, the error in fss/fs (given by the difference of the curves for these 
two quantities) can still be substantial. 

• There is a considerable dependence of the estimates on the unknown coupling gl. If 
we trust the large expansion, and assume \gi/ gql < 0.3, then we can draw further 
conclusions. In particular: 

— For small gl, the estimated errors in Bb and Bb, are small, less than roughly 0.1 
in magnitude. This is true for all values of g^ in our reasonable range. The error 
in the ratio BbJBb is smaller still, though of either sign. 

— The error in the ratio fss/fB is negative, and lies in range —0.3 to —0.1. As 
discussed above, this is because the quenched chiral logarithm decreases with 
rriss, due to the term proportional to 6. Taken literally, this would imply that the 
ratio in QCD is 10-30% larger than in QQCD. More conservatively, it suggests 
that the deviation of the ratio from unity, which in present QQCD simulations is 
ifss/ fs — 1) ~ 0.1 — 0.2, has a quenching error of 100% or more. 
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FIG. 6. Results for A{AX) for A = 1 GcV, a = 0, S = 0.15 and g ^ gg = 0.63. 

We show in Table |rV| our estimates of errors due to partial quenching. The results are 
for g"^ = QpQ = 0.4, but it is simple to scale to other values of g because the overall factors 
are the same for QCD and PQQCD. The results are fairly similar to those in Fig. ^ for 
g/ ^ 0. Most notably, the partial quenching error in fss/fB is negative and significant, just 
as for QQCD. The errors in 5-parameters are very small. 



V. CONCLUSIONS 



We have presented estimates for errors in lattice results for decay constants and B- 
parameters of heavy-light mesons due to quenching and partial quenching. Our estimates 
are reliable for the ratios Q/,_b, aside from higher order chiral corrections, or, in the case of 
PQQCD, corrections suppressed by m'jj/ml. Our estimates for the individual quantities, 
however, should be treated as approximate upper bounds. They suggest that, in order to 
extract fs and Bb, it is better to extrapolate from moderate quark masses than to attempt 
to work directly with very light valence quarks. As for the quenching errors themselves, 
those for 5-parameters are likely to be smaller than those in decay constants, with the 
latter possibly as small as a few percent, but more likely larger. Our most striking finding is 
that quenching reduces the ratio fss/fB, throwing doubt on the common assumption that 
this ratio is a fairly reliable prediction of QQCD. 
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FIG. 7. Results for A{AX) for A = 0.77 GeV, a = 0.7, 6 = 0.05 and g = gq = 0.63. 

Finally, we have found similar results for partially quenched QCD as for the fully 
quenched theory. This extends the observation of Booth P,p!0|, who found that there were 
significant differences between results for the unquenched theory with two flavors and QCD. 
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TABLE IV. Numerical results for A^^'^'^^ - A^'^^ , with = g^Q = 0.4. 



Qty. 


A = 


IGeV 


A = 0.77 GeV 






mjf = 0.15GeV2 


mjf = O.SGeV^ 


mjj- = 0.15GeV2 mjj- 


= O.SGeV^ 




-0.20 


-0.18 


-0.10 


-0.13 


fB 


-0.06 


0.05 


-0.04 


0.03 


IbJIb 


-0.14 


-0.24 


-0.06 


-0.16 


Bb. 


0.014 


0.010 


0.009 


0.008 


Bb 


-0.004 


-0.015 


-0.004 


-0.012 


BbJBb 


0.018 


0.025 


0.012 


0.020 
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